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The theory of geostrophic adjustment and its application to the dynamics of the high-latitude 
thermosphere have been discussed by ourselves (Larsen and Mikkelsen, JGR , 1983; Mikkelsen and 
Larsen, JGR , 1983) and Walterscheid and Boucher ( JAS , 1984) in previous papers based on a 
linearized treatment of the fluid dynamical equations. However, a linearized treatment is only 
valid for small Rossby numbers given by Ro = V/fL, where V is the wind speed, f is the local 
value of the Coriolis parameter, and L is a characteristic horizontal scale for the flow. For typical 
values in the auroral zone, the approximation is not reasonable for wind speeds greater than 25 
m/s or so. We have developed a shallow-water (one layer) model that includes the spherical 
geometry and full nonlinear dynamics in the momentum equations in order to isolate the effects of 
the nonlinearities on the adjustment process. A belt of accelerated winds between 60° and 70° 
latitude was used as the initial condition. We found that the adjustment process proceeds as ex- 
pected from the linear formulation, but that an asymmetry between the response for an eastward 
and westward flow results from the nonlinear curvature (centrifugal) terms. In general, the 
amplitude of an eastward flowing wind will be less after adjustment than a westward wind. For in- 
stance, if the initial wind velocity is 300 m/s, the linearized theory predicts a final wind speed of 
240 m/s, regardless of the flow direction. However, the nonlinear curvature terms modify the 
response and produce a final wind speed of only 200 m/s for an initial eastward wind and a final 
wind speed of almost 300 m/s for an initial westward flow direction. Also, less gravity wave 
energy is produced by the adjustment of the westward flow than by the adjustment of the 
eastward flow. The implications are that the response of the thermosphere should be significantly 
different on the dawn and dusk sides of the auroral oval. Larger flow velocities would be expected 
on the dusk side since the plasma will accelerate the flow in a westward direction in that sector. 
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Figure 1. The zonal wind (left), meridional wind (middle), and height of the fluid surface (right) at 
1 hr intervals. The initial condition at 0 hr is a 1 m/s eastward wind in the latitude band between 
60° to 70°. The height of the fluid is directly related to the pressure. All subsequent figures will be 
presented in the same format. 
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Figure 2. Analytic solution for the adjusted flow based on a linearized set of equations in a 
Cartesian coordinate system. The initial wind was 100 m/s, and the final wind is 80 m/s. The 
initial and final flow velocity will scale linearly together. Thus, an initial wind of 1 m/s will 
produce an adjusted wind of 0.8 m/s. This figure can be compared directly with the solutions 
shown in Figure 1. 
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Figure 3. An enlarged graph showing the adjusted state after 10 hr. The initial wind is eastward at 
1 m/s between 60° and 70 °N. The final zonal velocity agrees closely with the zonal velocity given 
by the linearized analytic solution. The meridional wind is negligible. The height of the free 
surface, which is directly related to the pressure, shows an asymmetry across the channel. The 
spherical geometry causes a focusing at the pole. As a result, the heights are higher on the 
southern side of the belt and lower on the northern side. 
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Figure 4. Adjusted state in the shallow-water model after 10 hr for an initial eastward wind of 100 
m/s. The final zonal wind is approximately 75 m/s which is less than the final value of 80 m/s 
based on the solution of the linearized set of equations. 
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Figure 5. The adjusted state for an initial 300 m/s eastward wind after 8 hr. The final zonal wind 
is close to 210 m/s which is nearly 30 m/s less than the predicted value based on a linear analysis. 
The oscillations which are evident in all three fields are due to gravity waves generated by the 
adjustment process. The waves propagate north and are reflected by the boundary condition at the 
pole. Therefore, the 300 m/s case is near the limit of validity of the numerical solution. 
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Figure 6. The adjusted state for an initial westward wind of 100 m/s is shown. The amplitude of 
the final zonal wind is approximately 90 m/s, or 10 m/s greater than expected based on a linear 
treatment. 
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Figure 7. Here the adjusted state for an initial westward wind of 300 m/s is shown. The final state 
also has a zonal wind of 300 m/s. The pressure gradient is very small compared to the cases 
shown in the previous figures, and there is almost no gravity wave energy generated, as can be 
seen in the smoothness and low values of the meridional wind. Thus, a significant asymmetry 
between the dusk and dawn sides of the auroral oval would be expected. We expect that on the 
dusk side, where the winds produced by the Lorentz acceleration are westward, less gravity wave 
energy will be produced by the adjustment process, and the pressure gradients will be smaller. 
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INITIAL WIND 
EASTWARD 

100 m/s 
200 m/s 
300 m/s 


LINEAR NONLINEAR 


80 m/s 80 m/s 

160 m/s 140 m/s 

240 m/s 200 m/s 


WESTWARD 


-100 m/s 
-200 m/s 
-300 m/s 


-80 m/s -90 m/s 

-160 m/s -200 m/s 

-240 m/s -300 m/s 


Figure 8. Table comparing the initial wind speed, the wind speed after adjustment in a linearized 
formulation, and the wind speed after adjustment when the nonlinear terms in the momentum 
equation are included. 
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